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Evolution of the fine-structure constant in the non-linear regime 
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We study the evolution of the fine-structure constant, a , induced by non-linear density perturba¬ 
tions in the context of the simplest class of quintessence models with a non-minimal coupling to the 
electromagnetic field, in which the two available free functions (potential and gauge kinetic function) 
are Taylor-expanded up to linear order. We show that the results obtained using the spherical infall 
model for an infinite wavelength inhomogeneity are inconsistent with the results of a local linearized 
gravity study and we argue in favour of the second approach. We also discuss recent claims that 
the value of a inside virialised regions could be significantly different from the background one on 
the basis of these findings. 


I. INTRODUCTION 

One of the deepest question of modern physics is 
whether or not there are fundamental scalar fields in 
nature. For example, they are a key ingredient in the 
standard model of particle physics (cf. the Higgs parti¬ 
cle, which is supposed to give mass to all other particles 
and make the theory gauge-invariant), but after several 
decades of accelerator experiments there is still no shred 
of experimental evidence for them. 

The early universe is a much better (not to mention 
cheaper) laboratory for fundamental physics. Observa¬ 
tions suggest that the recent universe is dominated by 
an energy component whose gravitational behaviour is 
quite similar to that of a cosmological constant (as first 
introduced by Einstein). This could of course be the 
right answer, but the observationally required value is so 
much smaller than what would be expected from parti¬ 
cle physics that a dynamical scalar field is arguably a 
more likely explanation. Now, the slow-roll of this field 
(which is mandatory so as to yield negative pressure) 
and the fact that it is presently dominating the universe 
imply (if the minimum of the potential vanishes) that 
the field vacuum expectation value today must be of or¬ 
der mpi , and that its excitations are very light, with 
to ~ H 0 ~ 10~ 33 eV. But a further consequence of this 
is seldom emphasized [jj: couplings of this field lead to 
observable long-range forces and time-dependence of the 
constant of nature (with corresponding violations of the 
Einstein Equivalence Principle). 

Measurements of various dimensionless couplings, such 
as the fine-structure constant a (which will be the focus 
of this paper) or the electron to proton mass ratio 3 
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are therefore unique tests of fundamental physics. Note 
that since the scalar field is effectively massless on solar 
system scales, it should in principle be easier to find new 
physics on astrophysical and cosmological scales. More¬ 
over, bounds on varying couplings restrict the evolution 
of the scalar field and provide constraints on dark en¬ 
ergy Si that, with new datasets becoming available 
in the near future will be complementary to (and indeed 
more powerful and constraining than) those obtained by 
traditional means. 

Let us now focus on the fine-structure constant, a , 
which among other things measures the strength of the 
electromagnetic interaction. The good news is then that, 
since the standard physics is changed in a number of key 
ways if there is a spacetime variation of a, there are many 
different ways in which measurements of a can be made. 
To name just a few, locally one can use atomic clocks || 
or the Oklo natural nuclear reactor ill On the other 
hand, on astrophysical and cosmological scales a lot of 
work has been done on measurements using quasar ab¬ 
sorption systems 8, 9. '.10j, Ilf and the cosmic microwave 
background [Til flil fill. I IT llfi|| . 

The bad news, however, is that these different mea¬ 
surements probe very different environments, and there¬ 
fore it is not trivial to compare and relate them. Sim¬ 
ply comparing at face value numbers obtained at differ¬ 
ent redshifts, for example, it is at the very least too 
naive, and in most cases manifestly incorrect. Indeed, 
detailed comparisons can often only be made in a model- 
dependent way, meaning that one has to specify a cosmo¬ 
logical model (crucial to define a clock in the universe, 
that is, a timescale) and/or a specific model for the evo¬ 
lution of a as a function of redshift. Simply assuming, 
for example, that alpha grows linearly with time (so that 
its time derivative is constant) is not satisfactory, as one 
can easily see that no sensible particle physics model will 
ever yield such a dependence for any significant redshift 
range. Last but not least, new methods are being devel¬ 
oped for measuring a using emission lines [fj], [h| [lf| as 
well as the electron to proton mass ratio (using absorp- 



2 


tion) m m m, so the issue of detailed comparisons 
between datasets will be even more important for the 
next generation of datasets. 

Here, we discuss one specific aspect of this issue. The 
scalar field responsible for the variation of a will (in any 
sensible particle physics model) couple to the matter sec¬ 
tor. Among other things, this implies that when, in 
the course of the cosmological evolution, inhomogeneities 
grow, become non-linear and decouple from the back¬ 
ground evolution, the same could happen to the local 
variations of ck. This has been previously studied in 
|23L I 24 I l2j using a simple spherical infall model for 
the evolution of infinite wavelength density perturbations 
and a particular generalization of the Bekenstein model 
for the evolution of a. (See also (2 tJ for a discus¬ 
sion of large scale variations of ck.) It was found that in 
the linear regime and in the matter era the variation of 
a would follow the density contrast. Moreover, it was 
also claimed that this approach was valid in the non¬ 
linear regime (meaning turnaround and collapse). Here 
we revisit these results, in particular questioning the ap¬ 
plicability of the spherical infall model. Note that in the 
particular case of the models of nail mill it is en¬ 
forced (purely by hand) that the evolution of a is driven 
by a coupling to charged non-relativistic matter alone 
and in that case the scalar field cannot also provide the 
dark energy. For that reason one should see these models 
as toy models that are useful for computational purposes. 
To some extent a similar comment applies to the models 
that we shall consider, though the reason here is purely 
their extreme simplicity. However, we do expect that at 
least at a qualitative level our results will be representa¬ 
tive of more realistic models. 

While this paper was being written up, ref. j2g] 
appeared. This provides a more detailed and 
mathematically-inclined analysis of local variations in 
physical ‘constants’, but does confirm our results. Our 
approach, while much simpler, has the advantage of mak¬ 
ing explicit the reasons why spatial variations have to be 
small, and why the use of the spherical collapse model of 
a infinite wavelength perturbation is inadequate. 

We will start in Sect. [©with a brief description of the 
models that we will be using. In Sect. IHIl we discuss the 
non-linear evolution of the fine structure constant using 
two different approaches. Finally, we describe and discuss 
our results in Sect. UY1 and present our conclusions in 
Sect. El Throughout this paper we shall use fundamental 
units with h = c = G = 1. 


II. THE LINEARIZED MODELS 


where 

£ 0 = icP 0 c> M 0 -I/( 0 ), ( 2 ) 

C+F = ~B f (4>)F im ,FI* v , (3) 

and £ 0 ther is the Lagrangian density of the other fields. 
We will make the simplifying assumption that both V 
and Bp are linear functions of 0 , namely 

1/(0) = V(0 O ) + ^ (0 ^ 0 O ) , (4) 

where dV/d<p is assumed to be a constant, the subscript 
‘ 0 ’ indicates that the variable local value is to be evalu¬ 
ated at the present time, and 

Bf{4>) = 1 - <fc(0 - 0o), (5) 

where k 2 = 8n and Equivalence principle tests require 
that iilU 

|C| < 5 x 10 -4 . ( 6 ) 

Since Bp = ao/a , this implies that 

— = 1 + (k (0 - 0o) , (7) 

CKO 

up to linear order. The equation of motion for the field 
0 is given by 

' < 8 > 

Therefore the evolution of the average value of ck with 

physical time is given approximately by 

a + 3Ha ~ {-£ip m + 6)«, (9) 

where we have defined F^F^/4 = — 7 Fp m , £1 = 
7f(C ^) 2 < 6 x 10 - 6 7_f and £2 = (k(dV /d0). On the 
other hand, linearizing eqn. © to obtain a local static 
solution in a slightly perturbed Minkowski space we ob¬ 
tain 

V 2 a 

~ £ 1 5pm- ( 10 ) 

CK 

In the following, for simplicity we shall assume that 
£2 = 0 and we will define £ = d. Although a gener¬ 
alization to models with dV/dcj) ^ 0 is straightforward, 
our main results will not be dependent on this particular 
assumption. 


We will consider the class of models where a neutral HI. NON-LINEAR EVOLUTION OF a 

scalar field is non-minimally coupled (via a gauge kinetic 

function Bp) to electromagnetism, namely We will consider the evolution of the value of ck due 

to the growth of a uniform spherical matter inhomogene- 
(1) ity with a final size ro using two different approaches. 


"b £(j)F T Mother ; 
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In what we shall refer to as the local approximation , we 
solve the Poisson equation m outside the spherical dis¬ 
tribution of mass, M , to determine the spatial variation 
of a as 

8a _ a(r) — a M 

a a 4nr ’ 

identifying a with a(r = oo). Given that M = 
47T(5p m r 3 /3 we have 

- = -W. (12) 

a 3 

Note that this calculation assumes that the spherical in¬ 
homogeneity has a radius significantly smaller than the 
horizon in order for eqn. G3 to be a valid approximation 
locally. 

A different calculation was introduced in [3 El] where 
the evolution of a uniform (that is, infinite wavelength ) 
perturbation in a was studied. The evolution of two ho¬ 
mogeneous and isotropic universes was considered with 
the initial conditions set at some initial time ti deep in 
the matter era. Here, we shall consider a similar calcu¬ 
lation where the background universe in a spatially flat 
Friedman-Robertson-Walker universe containing matter, 
radiation and a cosmological constant whose dynamics is 
described by 


a 


a 


= -Hi 



(13) 


with f Iki = 1 — — f l r i — Hai = 0 and taking the initial 

time ti to be deep into the radiation era. We also consider 
a perturbed closed universe with Hf ~ Hi, = i} ri , 
^l\i = Ha i, H mi = T and H^ = An m ,. In 

this paper, for simplicity we include the contribution of 
the energy density of the field <f> in the value of Ha thus 
neglecting the contribution of the kinetic contribution 
to energy density of the field </>, which is nevertheless 
constrained to be small. Also, since any variation of the 
fine structure constant from the epoch of nucleosynthesis 
onwards is expected to be very small PHI we neglect 
the minor contribution that such a variation has in the 
evolution of the baryon density (included in fi m o). 

In this approach the average values of a in the back¬ 
ground and perturbed universes are computed using eqn. 
@ so that 



=-) + 3H M (14, 


where 


Aa a — ai 
a di 


(15) 


and the subscripts b and oo indicate that a is calculated 
for the background and perturbed universes respectively. 



FIG. 1: Evolution of the scale factor, a, as a function of 
physical time, t, for the background and perturbed universes 
(solid and dashed lines respectively) 


By contrast, in the local approximation Aa/a can be 
computed as 


where 

5a 

a 



(16) 


^(H 2 n m ) p - (H 2 n m ) 

^ 8irH 2 


(i a p r 0 H) 2 . 


(17) 


We can already anticipate that these two different ap¬ 
proaches to estimate spatial variations of a will produce 
very different results since in the second approach there 
is no reference to the size of the fluctuation (an infinite 
wavelength approximation is considered). 


IV. RESULTS AND DISCUSSION 

The initial values of the various cosmological parame¬ 
ters were chosen in such a way that at the present time we 
have fl° m = 0.29, H Ao = 0.71 and H r0 = 8.4 x 1(T 5 jUJ 
for the background universe. Also, the value of Afl m ; 
is such that the collapse of the perturbed universe oc¬ 
curs near the present epoch. This is clearly seen in Fig. 
1 which plots the evolution of the scale factor, a as a 
function of the physical time, t, for the background and 
perturbed universes (solid and dashed lines respectively). 
Note that, according to the spherical collapse model, a 
uniform density perturbation will become infinite in a fi¬ 
nite time. Clearly this is not realistic (at least for local 
fluctuations) since in practice the collapse is stopped by 
virialisation. 

In Fig. 2 we compare the background evolution of 
Aa/a (dotted line) with that obtained using the two dif¬ 
ferent models described in the previous section. While 










4 



FIG. 2: Evolution of the fine structure constant, a, as a func¬ 
tion of physical time, t, in the vicinity of spherical distribution 
of mass including the spatial variations of a. calculated using 
the local approximation with ro = 2 Hq 1 (solid line) and the 
infinite wavelength approximation (dashed line). The dot¬ 
ted line represents the background evolution of a. Note that 
( < 6 x 1CP 6 7f is a very small number. 



FIG. 3: Evolution of (a — a/)/[a 6 — Qi), in various models 
as a function of physical time, t. The dashed line represents 
the infinite wavelength approximation while the solid lines 
use the local approximation with ro = 2 Hq 1 , ro = H// 1 and 
ro = Hq 1 /2 (top to bottom respectivelly). 


the evolution of A a/a obtained using the local approx¬ 
imation (dashed line) depends on the final size of the 
perturbation, ro, no reference to the perturbation size is 
made in the second case (solid line) where an uniform, 
infinite wavelength was considered. Hence, the latter 
model will produce a single result while the first method 
will produce a result which depends explicitly on the fi¬ 


nal perturbation size, ro- In Fig. 2 the results for the 
local approximation assume that ro = 2Hq . This is 
clearly an unrealistic value of ro which we use in order 
to get a sizable depart from the background A a/a (note 
that in order to be self-consistent the model requires that 
r <C .ff -1 ). However, the result can easily be rescaled for 
realistic values of ro using (IT71) (for example for a typical 
size of a galaxy cluster). 

In Fig. 3 we show the evolution of (a — ab)/{ab — 07 ), 
in various models as a function of physical time, t. The 
dashed line represents the infinite wavelength approxi¬ 
mation while the solid lines use the local approximation 
with ro = ZHq 1 , ro = Hq 1 and ro = Hq 1 /2 (top to 
bottom respectivelly). We clearly see that as we move 
towards smaller (more realistic) values of ro the local ap¬ 
proximation will give negligible variations of a. 

It is interesting to discuss the evolution of 5a/a in 
the linear regime during the matter and radiation eras. 
Note that in the linear regime the physical radius of the 
spherical overdensity grows as 

r oc a , (18) 

and therefore 


M oc 6p m a 3 = 6mPma 3 


(19) 


which is approximately constant during the radiation era 
and grows proportionally to a during the matter era. 
Hence deep into the radiation era 


5a , 

— oc a , 

(20) 

a 

while deep into the matter era 


5a 

— ~ constant. 

(21) 


a 

It is also important to point out that the amplitude of 
the a variation (for a given 5p m ) depends quadratically 
on the radius of the spherical overdensity (Aa/a oc r 2 ). 
Consequently, for ro <C Hq 1 the spatial variations pre¬ 
dicted by the local approximation will be very small. The 
background evolution of A a/a in linearized Bekenstein 
models has been previously studied in [§2] and it is well 
known that A a/ a —> 0 deep into the radiation era. 

This is confirmed by the results plotted in Fig. 2 which 
clearly confirm the above discussion. In both approxima¬ 
tions 


Aa 

a 



( 22 ) 


deep into the radiation era. On the other hand 


~ constant, (23) 

b 

deep into the matter era. Also, as expected the evolution 
of Aa/a in the infinite wavelength approximation only 
departs significantly from the background one near the 
present time (when the collapse of the perturbed universe 
occurs). 
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V. CONCLUSIONS 

We have presented and contrasted two different ap¬ 
proaches for the determination of the non-linear evolu¬ 
tion of a. While in the local approximation we assume 
that the size of the perturbation which gives rise to spa¬ 
tial variation of a is much smaller than the Hubble radius, 
if -1 , in the infinite wavelength approximation the oppo¬ 
site is required for self-consistency. Hence, it is not sur¬ 
prising that the results obtained using both models are so 
different. The crucial question is which one provides the 
right answer when applied to cosmology. Here, we argue 
that local approximation provides the right answer since 
the non-linear effects are only expected to be important 
in this context on scales much smaller than H^ 1 (scales 
smaller than the a typical galaxy cluster size). On such 


small scales, we therefore predict that the spatial vari¬ 
ations of a generated in the simplest models should be 
too small to be of any cosmological interest. These re¬ 
sults confirm that it is difficult to g enerate significant 
large-scale spatial variations of a [27], [33, j34| even when 
we account for the evolution of the fine structure constant 
in the non-linear regime. 
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